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Abstract. The gravitational wave (GW) signature of a binary black hole (BBH) coalescence 
is characterized by rapid frequency evolution in the late inspiral and merger phases. For a 
system with total mass larger than 100 Mq, ground based GW detectors are sensitive to the 
merger phase, and the in-band whitened waveform is a short-duration transient lasting about 
10-30 ms. For a symmetric mass system with total mass between 10 and 100 Mq, the detector 
is sensitive instead to the inspiral phase and the in-band signal has a longer duration, between 
30 ms -3 s. Omega is a search algorithm for GW bursts that, with the assumption of locally 
stationary frequency evolution, uses sine-Gaussian wavelets as a template bank to decompose 
interferometer strain data. The local stationarity of sine-Gaussians induces a performance loss 
for the detection of lower mass BBH signatures, due to the mismatch between template and 
signal. We present the performance of a modified version of the Omega algorithm, Chirplet 
Omega, which allows a linear variation of frequency, to target BBH coalescences. The use 
of Chirplet-like templates enhances the measured signal-to-noise ratio due to less mismatch 
between template and data, and increases the detectability of lower mass BBH coalescences. 
We present the results of a performance study of Chirplet Omega in colored Gaussian noise at 
initial LIGO sensitivity. 



1. Motivations 

Gravitational wave burst searches [1] are designed to search for short duration transients. 
Burst search algorithms typically look for excess signal energy in the time-frequency plane 
after filtering the GW strain data with a particular template bank. Current burst algorithms, 
such as Omega [2] and Coherent Waveburst [3] use sine-Gaussians and a collection of Meyer 
wavelet bases respectively as templates. Matched filtering searches for the coalescence of compact 
binary systems with total mass 2 — 25Mq [1], 25 — WOMq [5] and > IOOMq [6] have also been 
performed. These matched filtering algorithms use signal models such as the Post Newtonian 
approximation [7j, the EOBNR analytical template family [8] and quasinormal ringdown modes, 
respectively. 

Inspiralling binary black holes are characterized by chirping GW signals, where the frequency 
changes rapidly in the later inspiral and the merger phases. The duration of the signal in the 
most sensitive band of a ground based detector is dependent on the mass of the binary. For 
a system of total mass larger than 100 Mq, a ground based GW detector is sensitive to the 
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(a) A sample sine-Guassian. 
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Figure 1: Sine-Gaussian and Chirplet waveform: (left) Example of a Sine-Gaussian 
waveform, (right) Example of a Chirplet waveform. 



merger phase, and the in-band whitened waveform is a short-duration transient lasting 10-30 
ms. For a symmetric binary system with a total mass between 10 and 100 Mq, the detector is 
sensitive instead to the inspiral phase and the in-band signal has a longer duration between 30 
ms - 3 s. A chirplet based burst search has been proposed in Ref. [9], which combines the chirp 
characteristic of an inspiral with the existing Omega burst search. In this paper we explore 
further the use of the chirplet burst search and benchmark the detection performance for BBH 
coalescences, compared to the standard Omega burst search at a constant false alarm rate. 

2. Chirplet Definition 

The sine-Gaussian template bank used in the standard Omega burst search, are defined in the 
time domain as: 



V'(t) = Aexp 



(r-t)^ exp(2^i [f{T-t)]), 



Chirplets are a generalization of sine-Gaussians, defined in time domain as: 

(2vr/)2 



Aexp (-^^(^ - tf) exp (2vri [/(r - t) + d/2 (r - tf]) 



(1) 



(2) 



where A = (87r/2/Q2)i/4 

is a normalization factor ensuring that / I'lpl'^dT = 1; t and / are the 
central time and central frequency, respectively and Q is the dimensionless quality factor. See 
Fig. [Hfor an example of sine-Gaussian and a chirplet. 

The difi^erence between a chirplet and a sine-Gaussian waveform is the chirp rate, an additional 
term in the phase, denoted as d, that linearly controls the slope of the chirplet frequency evolution 
as /(r) = / + d{T — t). When d = 0, a chirplet becomes the standard sine-Gaussians. Chirplets 
are thus described with a four-dimensional parameter space of {t, f,Q,d}. 

The chirplet transform T is obtained by correlating the data with chirplets. In the frequency 
domain, it reads: 



T[x-{tJ,Q,d}] 



XiC)^f*i^;{t,f,Q,d})dC 



(3) 



where X{S^) and ^'(^; {t, f,Q,d}) denote the Fourier transforms of the (whitened) data stream 
x(r) and chirplet 'iP{t) with the parameters {t,f,Q,d} respectively. 
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(a) Eventgram in sine-Gaussian bank. 
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(b) Eventgram in chirplet banlt. 



Figure 2: Eventgrams for sine-Gaussian and Chirplet banks: The eventgrams are plotted 
for inspiralhng black- hole binary (with mi = 13Mq and m2 = 20Mq) signal in simulated 
Gaussian LIGO noise. See the difference between the constant frequency tiles for the sine- 
Gaussians and the linear frequency variation for chirplets. An increase of SNR by ~ 20% in the 
peak pixel (the most significant time-frequency pixel) is to be noted. 



3. Performances of Chirpletized Omega 

3.1. Searching for binary coalescence with chirplets 

The gravitational wave spectrum from the coalescence of a binary system is at lower frequency for 
system with larger total mass. In particular, the frequency associated with the inner-most stable 
circular orbit (ISCO), which corresponds to the transition between inspiral and merger phase 
of the coalescence, is below 80 Hz for systems with total mass M > 100-Mq. With this feature, 
the chirp phase of the waveform with a spectral content at frequencies below ISCO is outside 
the detector sensitive band and thus does not contribute significantly to the total detected SNR 
of the signal. Sine-Gaussian waveforms are a reasonable approximation for the few waveform 
cycles associated with the merger and ringdown portions of the coalescence waveform. For a 
symmetric system with lower mass {M < 30Mq), the coalescence time of BBH in a ground 
based detector's sensitive band is > 400 ms. There is loss of measured signal-to-noise ratio 
induced by the mismatch in the waveforms; between the chirp nature of the GW evolution and 
the constant central frequency sine-Gaussian template bank. The linear frequency variation of 
chirplets provide better match to the chirping signal of BBH in this mass range. 

3.2. Increased signal-to-noise ratio 

It has been shown [9] that Chirpletized Omega measures a higher matched filtering signal-to- 
noise-ratio (SNR) for BBH signal < IOO-Mq. The SNR improvement is (~ 20%) for the mass 
range M < 60Mq and the improvement may go up-to ~ 40% for M < 20AfQ. When compared 
to actual matched filtering SNR of signals, chirplet recovers the full SNR for M < 2O-M0. 

As an illustration, we show an eventgram (Fig. [2j). An eventgram shows the significant time- 
frequency pixels for different Q values. In this particular eventgram a BBH signal is added to 
Gaussian noise with the spectral characteristics of initial LIGO design. The signal is a non- 
spinning BBH system with mass components mi = 13Mq and m2 = 2O-M0. In this example, 
chirplets with a positive slope are chosen to sine-Gaussian with constant frequency and the 
correlation of the most significative chirplet is, ~ 20% larger than the most significative sine- 
Gaussian. 
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(b) Efficiency curves for set 35-80 Mg 



Figure 3: Efficiency curve: Efficiency sigmoids are plotted with a FAR of 10~^ Hz. Higher 
detection efficiency in a chirpiet search is observed. The improvement in chirpiet is more 
pronounced for the set with mass 4-35 Mq. 



3.3. Increased detectability 

The detectability of BBH signals is limited by a background of noise fluctuations that trigger 
the search algorithms. Noise fluctuations can be reduced by eliminating non-GW triggers from 
the search with the application of a suitable coincident or coherent strategy for a network of 
GW detectors. The background in a chirpiet search is in general different from the typical 
background in the standard Omega burst search. An additional 5% increase was reported [9] in 
the measured single detector SNR for background triggers in a chirpiet burst search compared to 
standard Omega. Previous analysis was conducted with different false alarm rates. In order 
to ensure an unbiased comparison, this paper addresses this issue and compares the detectability 
in the two search algorithms at the same false alarm rate (FAR) . In the current study we added 
BBH signals to 2 months of simulated colored Gaussian noise mimicking the instrumental noise 
of the two Hanford LIGO detectors (HI and H2), and the Livingston LIGO detector (LI). We 
split the range of BBH masses in two sets, total mass 4-35 Mq and total mass 35-80 Mq. The 
distribution of mass was uniform in the component masses of the two black holes. We estimated 
the network's background by time-sliding the triggers from LI detector with the coincident 
triggers from HI and H2. Similarly foreground triggers were collected by coincidence of triggers 
from HI, H2 and LI. A time window of 10 ms was used to find the coincident triggers. The 
times at which the signals were added to noise are known. We used those times to identify 
significant events within 10 ms in the foreground triggers. We ranked the triggers from both 
the background and the simulated signals using the the geometric mean of the single detector 
signal-to-noise ratios as ranking statistics. 

1 

PGM = {PHl X PH2 X /OLl)S (4) 

where pdetector is the single detector SNR. 

After all the triggers are ranked, the significance of each trigger was evaluated by computing 
its probability to belong to background (FAR). The performance comparison between chirpiet 
and standard Omega searches was done with a common threshold on the FAR. 



Efficiency curve — We compared the performance of the sine-Gaussian and the Chirplet 
algorithms by estimating the detection efficiency as a function of the effective distance of the 
BBH from the detector. We fit a sigmoid to the observed number of signals as a function of the 
distance, at a fixed false alarm rate of 10~^ Hz. The characteristic distance parameters, known 
as the 50% and 90% efficiency distances, which are the points on the sigmoid where 50% and 
90% of the signals are recovered respectively, are both improved for chirplet compared to Omega 
as seen in Fig. [3J The numbers obtained from the sigmoid fit are shown below in Table 1. The 
overall improvement of > 50% in the distance parameters, is higher for the lower mass set (4-35 
Mq) compared to, up-to 40% for the mass set with total mass ranging 35-80 Mq. 



Set 


Sine-Gaussians 
50% (Mpc) 


Sine-Gaussians 
90% (Mpc) 


Chirplets 
50% (Mpc) 


Chirplets 
90% (Mpc) 


4-35 Mq 


56 ± 4 


19 ± 1 


84 ± 2 


35 ± 2 


35-80 Mq 


139 ± 4 


66 ± 2 


166 ± 4 


75 ± 3 



Table 1: 50% and 90% efficiency distance at a FAR of 10"^ Hz. 



Sensitivity distance — While efficiency sigmoids are good measures for the detection 
performance, as shown in Fig. [3l they don't provide details of how detectability changes across 
mass bins. For this reason, we compared the fiducial distance the detector network is sensitive 
at a false alarm rate of 10~^ Hz, or the sensitivity distance, as the function of the BBH mass 
parameters. Fig. [Ja] and Fig. I4bl show sensitivity distance for chirplet and standard Omega 
search in the mass parameter space for the set with total mass 4-35 Mq . Fig. He] shows the ratio 
of the sensitivity distance measured by chirplet to the sensitivity distance measured by standard 
Omega search. For this set, chirplet search improves detectability more than 30% over standard 
Omega search in all mass bins as seen in the ratio plot of Fig. HcJ Similarly, Fig. I4dl and Fig. Hel 
show detectability for set with total mass 35-80 Mq. In this set the detectability in chirplet 
search also improves up-to 30% over standard Omega search for most of the mass bins and is 
shown in Fig. |4fl 

The above comparisons of efficiency sigmoids and sensitivity distances at a false alarm rate 
of 10~* Hz point very well to higher detectability of BBH signal in a chirplet burst search over 
standard Omega burst search, specially in the lower mass region (total mass < 35 Mq). This is 
in accordance with the previous result of increased SNR [S] . 

4. Future plans 

In this study, we have demonstrated the advantage of a chirplet burst search over the 
standard sine-Gaussian burst search using coincidence method as the multi-detector network 
strategy. There is room for additional modifications in the chirplet-based analysis, with the 
implementation of a coherent detection statistics and novel clustering strategies, which are 
currently under study. There have been studies \10\ [TT] based on matrix factorization techniques 
to construct orthonormal bases that span the space of high-mass merger waveforms. Using 
chirplets as templates, is a similar but distinct approach, with different advantages. In particular, 
a chirplet based search is not tied to the availability of an explicit representation of the family 
of waveforms to construct the basis, which gives it greater simplicity and robustness in the face 
of possible errors in the numerical relativity simulations, although, in exchange, it suffers some 
possible SNR loss. 
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Figure 4: Sensitivity distance: Plots of the fiducial distance that can be observed at a FAR 
of 10~^ Hz for both chirplet and sine-Gaussian in the mass parameter space are shown in Fig. 
3al I4bl I4dl and Hel The ratio plots of the detectable distance are shown in Fig. He] and More 
than 30% of improvement in the sensitivity distance for the set 4-35 Mq in all the mass bins 
and up-to 30% improvement in the sensitivity distance for most of the mass bins in set 35-80 
Mq should be noted. 



In addition to an improved detectability, the chirplet-based analysis offers additional 
information on the frequency evolution of the observed signal, which will be useful in the a 
posteriori interpretation of an event. 

The work presented in this paper does not address the question of the spin of the binary 
black hole system. Spins affect the detectability profoundly [12]. We are planning to extend our 
current studies to express the detectability measured in a chirplet burst search to spin parameter 
space. 
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